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Abstract 

We draw a possible scenario for the observation of massive long-lived charged particles at the 
LHC detector ATLAS. The required flexibility of the detector triggers and of the identification 
and reconstruction systems are discussed. As an example, we focus on the measurement of the 
mass and lifetime of long-lived charged sleptons predicted in the framework of supersymmetric 
models with gauge-mediated supersymmetry (SUSY) breaking. In this case, the next-to-lightest 
SUSY particle can be the light scalar partner of the tau lepton (fi), possibly decaying slowly 
into a gravitino. A wide region of the SUSY parameters space was explored. The accessible 
range and precision on the measurement of the SUSY breaking scale parameter y/~F achievable 
with a counting method are assessed. 
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1 Introduction 



The Large Hadron Collider (LHC) at the European Laboratory for Particle Physics 
(CERN) in Geneva has an enormous and unique potential of discovering new particles 
beyond the Standard Model (SM) over a large mass range around the TeV scale. 
In this article, we consider an interesting class of exotic particles that are characterised 
by a few model independent common features: 

a) they have large masses and are produced with (3 < 1 in a non negligible fraction of 

events at the LHC; 

b) they are stable or quasi-stable, i.e. their lifetime is larger than 10~ 7 s, so that they 

decay far from the collision point, possibly outside the detector; 

c) their electric charge is either an integer or a fractional multiple of the proton charge. 

Actually, quite a few theories beyond the SM foresee charged particles, both strongly 
and electroweakly interacting, with large masses and, although many exotic particles are 
unstable, often the lowest (or the next-to-lowest) lying state may be stable or quasi- 
stable [|l],|2|. In the following, we call these states Massive Semi-stable Exotic Particles 
(MSEP's). 

In the first part of this paper, we will briefly recall the difficulties which a large apparatus 
like ATLAS could meet while trying to detect MSEP's. In the second part, we show a 
concrete example of the measurements (involving detection, identification and track re- 
construction) of MSEP's coming from supersymmetric models with gauge-mediated 
supersymmetry breaking (GMSB) In this scheme, the supersymmetry (SUSY) || 

breaking occurs at relatively low energy scales, and it is mediated mainly by gauge interac- 
tions. The automatic suppression of the SUSY contributions to flavour- changing neutral 
currents and CP-violating processes is naturally fulfilled. Furthermore, in the simplest 
versions of GMSB, the Minimal Supersymmetric Standard Model (MSSM) spectrum and 
other observables depend on just five parameters, usually chosen to be ||7|-[T0|: the over- 
all messenger scale M mess , the so-called messenger index N mess , the universal soft SUSY 
breaking scale felt by the low-energy sector A, the ratio of the vacuum expectation values 
of the two Higgs doublets tan/3, and sign(/x), which is the ambiguity left for the SUSY 
higgsino mass after imposing the conditions for a correct electroweak symmetry breaking. 

In the following, we will consider GMSB models where the role of MSEP is taken by 
the stau (or all charged sleptons), which is the next-to- lightest SUSY particle (NLSP), 
decaying into gravitinos. These scenarios are very promising at the LHC, providing signa- 
tures of semi-stable charged tracks coming from massive sleptons, therefore with f3 often 
significantly smaller than 1. In particular, we perform our simulations at the ATLAS 
muon detector, whose excellent time resolution Jl l] allows a precision measurement of 



the slepton time of flight, and hence of the slepton velocity. We show that the event 
can be recognised with high efficiency by the calorimetric trigger requiring the usual 
SUSY signature of P™ lss +jets, where the hadronic jets, originated by the decay of heavy 
squarks/gluinos along with MSEP's, must have the transverse momentum Pt > 50 GeV 
and the transverse momentum imbalance P™ 1SS is calculated only from the energy deposi- 
tion in the hadronic calorimeter. It is remarkable that, in this framework, the event rate 
is very high and the background contamination due to SM channels can be reduced to a 
completely negligible level. Then, we conclude stressing the fact that the measurements of 
mass (using the MSEP momentum obtained by the muon system) and lifetime (obtained 
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by a counting method) are sufficient to determine the SUSY breaking scale \/~F at a level 
of precision of a few 10%'s. 

Although this topic did not receive much attention up to now, as we will discuss, there 
is a sound theoretical basis to emphasise the physics and experimental search for MSEP's. 



2 The Experimental Challenge 

The experimental signatures of heavy long-lived charged particles at a hadron collider 
have been studied both in the framework of GMSB and in more general scenarios [|T^-[T^|. 
The two main observables one can use to separate these particles from muons are the 
high specific ionisation and the time of flight in the detector. As a matter of fact, a 
MSEP behaves like a massive muon with a velocity sometimes considerably lower than 
c. The energy loss of a MSEP in matter has been carefully studied in Ref. JT5J (see also 
Ref. where it was shown that, for energy loss for ionisation when (3 > ~ and for 
particles with M >100 GeV/c 2 , the range in iron exceeds several metres and increases 
with M. Moreover, for strongly interacting MSEP's the penetration length in matter 
does not change significantly for f3 < 0.7 and M =100 GeV/c 2 |TJ|. At the same time, a 



remarkable property of MSEP's is that, because of their low (3 value, they can produce 
anomalous energy loss for ionisation in matter that can be used as a distinctive feature 
for their identification. 

Nevertheless, there are some drawbacks associated to low-velocity MSEP's. In fact, 
the detector synchronisation is tuned for particles travelling across the apparatus with 
the velocity of light, hence, due to the large geometrical dimensions of the typical LHC 
apparatus, all the data acquisition system should be carefully tuned in order to guarantee 
the trigger and the event reconstruction of a MSEP. The trigger system of the ATLAS 
experiment is described in detail in Ref. (see Figs. [I] and |2] for schematic views of the 
ATLAS apparatus). Three levels of trigger are envisaged, which will reduce the initial 
input rate of 40 MHz for the first-level trigger progressively down to a rate of ~ 100 
Hz. This is the maximum input rate of the data storage system, and corresponds to an 
acquisition of ~ 100 MB/s. The first level comes exclusively from the hardware, and the 
information from the muon detectors and from the calorimeters are treated separately. 
The second level refines the first level by connecting the information from different detec- 
tors. Finally, the third level (also called event filter) applies the full off-line reconstruction 
algorithms to the data. 

The track identification of a MSEP escaping the hadronic calorimeter in the high 
muon region is the natural task of the ATLAS specialised structure of muon precision 
chambers, composed by MDT's (Monitored Drift Tubes), which are dedicated to mea- 
sure the trajectory of penetrating particles and to determine, with a good precision, the 
position of a charged particle. 

We recall that in the central region an air-core toroidal superconducting magnet gener- 
ates a magnetic field of about 0.5 T to allow the measurement of the particle momentum. 
The precision chambers are composed by two multilayers, and each of them is composed 
of three (four in the inner station) drift tube layers. A typical track goes through about 
twenty MDT's, before escaping from the apparatus. A single MDT is a cylindrical drift 
tube of 3 cm diameter. It is basically filled with Argon + CO2 at 3 atmosphere absolute 
pressure. It has a typical spatial resolution a = 80 fim and a maximum drift time of 700 



ns. The resulting spatial resolution for the single particle is about 30 /zm [[17]]. Once an 



event has been accepted by the first level trigger, all the data coming from the MDT's for 
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Figure 1: Isometric view of the ATLAS detector with the description of the different 
functional parts. 



the following 700 ns are labelled with the same bunch crossing number and are extracted 
from the pipe-line memories. The track reconstruction for a MSEP will be done by taking 
into account in the hit analysis the delay along the track due to (3 < 1. 

In our case, we will assume to trigger the MSEP's event by the hadronic calorimeter 
which can be easily activated by the multiple high Pt jets produced by the decay of 
squarks/gluinos with masses larger than 500 GeV together with MSEP's with masses 
of the order of 100 GeV. The classical SUSY signature of P™ lss +jets, where P™ 1SS is 
calculated only from the energy deposit in the calorimeter, neglecting the NLSP's and the 
muons, is in fact enough to pass the requests of the first level P™ 1SS trigger of a jet with 
Pt > 50 GeV and P™ 1SS > 50 GeV. A drawback of this purely calorimetric approach is 
the fact that processes with low hadronic activity, such as direct slepton production and 
direct electroweak gaugino production are not selected. 

The measurements of the time of flight for MSEP's are made possible by the timing 
precision ( < 1 ns) and the size of the ATLAS muon spectrometer. 

In the barrel part of the detector < 1), the precision muon system consists of three 
multilayers of precision drift tubes immersed in a toroidal air-core magnetic field. The 
three measuring stations are located at distances of approximately 5, 7.5 and 10 m from 
the interaction point. A particle crossing a drift chamber ionises the chamber gas along 
its path, and the electrons produced by the ionisation drift to the anode wire under the 
influence of an electric field. The particle position is calculated from the measurement 
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Figure 2: Side view of the ATLAS detector with the description of the different functional 
parts. 



of the drift time of the ionisation electrons to the anode wire. In order to perform this 
calculation, a starting time t for counting the drift time is needed, corresponding to 
the time of flight of the particle from the production point to the measuring station. 
For a particle travelling approximately at the speed of light, as a muon, the to's for the 
measuring stations are parameters of the detector geometry and of the response of the 
front-end electronics ||18| . For a heavy particle, the to is a free parameter, function of the 
(3 (= v/c) of the particle. It was demonstrated with a full simulation of the ATLAS muon 
detector [TS[] that the (3 of a particle can be measured by adjusting the to for each station 
in such a way to minimise the \ 2 °f the reconstructed muon track. 

The particle f3 can be measured with a resolution approximately parameterised as 
a(f3)/(3 2 = 0.028. The resolution on the transverse momentum measurement for heavy 
particles is found to be comparable to the one expected for muons. 

Besides the strong signature of the delay needed to reconstruct the track, a MSEP is 
identified by the energy loss for ionisation that can be much higher than that of a muon 
at minimum. This signature is a powerful physical quantity that gives a clear MSEP's 
identification. 

[12| , already at (3 



As discussed in Ref. 



0.6 the energy loss is around twice the 
minimum, and of course it increases rapidly for a lower value of (3. A complete discussion 
of the ionisation measurement in the ATLAS detector is a very complicated problem 
which can not be addressed in this work. Here, we want to illustrate only in a preliminary 
way the possible role of the hadronic calorimeter and of the MDT's. We believe that 
the calorimeter can bring a significant contribution to the identification of MSEP's. The 
ATLAS hadron calorimeter can measure the muon energy deposit with an accuracy of 
about 25%, in the last compartment, and ~ 20%, using the full depth. The contamination 
induced by the minimum bias event pile-up is expected to be of the order of 1% or lower 
at full luminosity. 

In the case of the MDT detector, the measurement of total charge based on a single 
MDT appears poor because of the background contamination at full LHC luminosity. 
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Notwithstanding, measurement of charge performed by a statistical procedure involving 
the complete set of ~ 20 tubes fired by the particle can provide a significantly improved 
measurement. It must be considered that the MDT's front-end electronics is optimised 
to provide a charge integration only for a small fraction of the total drift time (~ 30 ns). 
However, even with this limitation, it could be possible that again the statistical combi- 
nation of all the tubes belonging to the track can provide a good ionisation measurement. 
Of course, this brief discussion does not exhaust all the possibilities offered by the 



ATLAS apparatus. For example, it was shown in Ref. JTTJ] that the ionisation energy loss 
measurement in the Transition Radiation Tracker can be used to achieve ir/K separation 
and, in some cases, the electromagnetic calorimeter can provide the charge measurements. 



3 Supersymmetric Scenario 

The phenomenology of GMSB (and more in general of any theory with low-energy SUSY 
breaking) is characterised by the presence of a very light gravitino G 



"T.3/2 = m G 



V3M'r 



where \JF is the fundamental scale of SUSY breaking. If we assume the typical value 
y/F = 100 TeV, and M' P = 2.44 x 10 18 GeV for the reduced Planck mass, we get m 3 / 2 = 
2.37 eV. Hence, G is always the lightest SUSY particle (LSP) in these theories. If im- 
parity is assumed to be conserved, any produced MSSM particle will finally decay into 
the gravitino. Depending on \/F, the interactions of the gravitino, although much weaker 
than gauge and Yukawa interactions, can still be strong enough to be of relevance for 
collider physics. As a result, in most cases the last step of any SUSY decay chain is the 
decay of the next-to-lightest SUSY particle (NLSP), which can occur outside or inside 
a typical detector (even close to the interaction point). The pattern of the resulting 
spectacular signatures is determined by the identity of the NLSP and its lifetime before 
decaying into G, 



ctnlsp 1 / V F \ f m NLSP \ 5 
cm 100S VlOOTeV/ V 100 GeV J ' KJ 

where B is a number of order unity depending on the nature of the NLSP. The identity of 
the NLSP [or, to be more precise, the identity of the sparticle(s) having a large branching 
ratio (BR) for decaying into the gravitino and the relevant SM partner] determines four 
main scenarios giving rise to qualitatively different phenomenology: a) the Neutralino- 
NLSP Scenario; b) the Stau-NLSP Scenario; c) the Slepton co-NLSP Scenario; d) the 
Neutralino-Stau co-NLSP Scenario . In this paper, we will be concerned with Scenarios 
b) and c) only. 

The two main parameters affecting the experimental measurement at the LHC of the 
slepton NLSP properties are the slepton mass and momentum distribution. Indeed, at a 
hadron collider most of the NLSP's come from squark and gluino production, followed by 
cascade decays. Thus, the momentum distribution is in general a function of the whole 
MSSM spectrum. However, one can approximately assume that most of the information 
on the NLSP momentum distribution is provided by the squark mass scale rrtq only (in the 
stau NLSP scenario or slepton co-NLSP scenarios of GMSB, one generally finds rrig > rriq). 
To perform detailed simulations, we select a representative set of GMSB models generated 
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528420 


7437 
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BKGD: Iv+Jets 


9.6 


9.6 


9.5 


9.5 


O A 

Z.4 


O A 

z.4 


1.0 


1.8 


BKGD: Z+Jets 


6.8 


6.9 


6.9 


7.1 


11.1 


11.0 


5.9 


8.7 


BKGD: it 


5.3 


5.3 


5.3 


5.6 


6.2 


6.5 


3.4 


4.9 


BKGD: QCD 


8.0 


8.0 


8.0 


7.4 


3.1 


3.1 


0.5 


2.0 


BKGD: Total 


29.7 


29.9 


29.9 


29.9 


22.8 


23.0 


10.8 


17.4 



Table 1: In rows 2-5, the input parameters of the eight sample GMSB models chosen for 
our study are reported. In rows 6-10, the main features of the models (£ = eR, {Ir, fi) are 
shown. The average squark mass is indicated by rriq and rrig is the gluino mass. In the last 
six rows, the comparison among the expected number of events after the cuts described 
in Sect. 3.1 and the major background sources are reported. The assumed integrated 
luminosity is 30 fb , corresponding to a three-year low- luminosity run at the LHC. In 
the table, BKGD is a shorthand for "background" , and M mess and N mess are defined in 
the text (see Refs. 



by SUSYFIRE [|2T|. We limit ourselves to models with m-NLSP > 100 GeV, that cannot be 
excluded by direct searches at LEP/Tevatron, and rrtq < 2 TeV, in order to yield an 
adequate event statistics after a three-year low-luminosity run (corresponding to 30 fb -1 ) 
at the LHC. Within these ranges, we choose eight extreme points (four in the stau NLSP 
scenario and four in the slepton co-NLSP scenario) allowed by GMSB in the (ttinlsp, mq) 



plane, in order to cover the various possibilities. More details can be found in Refs. [p2| , |23[| . 
In Tab. [I] (first five rows), we list the input GMSB parameters we used to generate these 
eight points. 

3.1 Event Selection and Slepton Mass Measurement 

In order to select a clean sample of sleptons, we applied the following requirements: 

• at least a hadronic jet with Pt > 50 GeV and a calorimetric 
^miss > 5Q Q e y (trigger requirement); 

• at least one slepton candidate satisfying the following cuts: 

— |?7| <2.4 to ensure that the particle is in the acceptance of the muon trigger 
chamber, and therefore both coordinates can be measured; 
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— Aneas < 0.91, where /3 meas is the j3 of the particle measured with the time of 
flight in the precision chambers; 

— The Pt of the slepton candidate, after the energy loss in the calorimeters 
has been taken into account, must be larger than 10 GeV, to ensure that the 
particle traverses all the muon stations. 

a cut m c fj > 800 GeV, where m e g is the total invariant mass of the event constructed 
starting from the transverse momentum of the high Pt jets and muons (or muon-like 
particles) 



Considering an integrated luminosity of 30 fb _1 , a number of events ranging from a few 
hundred, for the models with 2 TeV squark-mass scale, to a few hundred thousand, for a 
500 GeV mass scale, survives these cuts. These events can be used to measure the NLSP 
properties. For the sake of comparison, we recall that about 1500 events for a light Higgs 



boson production are expected in the channel H — > 77 in the same environment [JTT . 

In order to perform the mass measurement, the particle momentum is needed. The 
precision chambers only provide a measurement of the momentum components transverse 
to the beam axis, so a measurement of the slepton pseudorapidity is necessary. This can 
be performed either by a match with a track in the inner detector, or using the information 
from the muon trigger chambers. The first option requires a detailed study of the matching 
procedure between detectors. This study was performed for muons in Ref. |IT |, but the 



results can not be transferred automatically to the case of heavy particles for which the 
effect of multiple scattering in crossing the calorimetric system is much more severe. 

In the case of the trigger chambers, a limited time window around the beam crossing is 
read out, restricting the (3 range for which the momentum can be measured. We therefore 
evaluated the statistical precision achievable for the eight example models in two different 
P intervals: 0.6 < (3 < 0.91 and 0.8 < f3 < 0.91. 



Many more technical details about this analysis can be found in Refs. [22, 23 . 



3.2 Slepton Lifetime Measurement and \J~F 

The measurement of the slepton NLSP lifetime was performed by exploiting the fact that 
a couple of NLSP's is produced in each event. We adopted a method similar to the one 



used in Ref. |T0[ for the high energy e + e~ collider in the neutralino NLSP case. 

This method consists in selecting the Nx events where a slepton is detected through 
the time-of-flight measurement and then counting how many times {N 2 ) a second slepton 
is observed. This information is sufficient to determine the lifetime. Although in principle 
very simple, in practice this method requires an excellent control of all possible sources 
of inefficiency for detecting the second slepton. 

The first set of Ni events is defined with the additional requirement that, for a given 
value of the slepton lifetime, at least one of the produced sleptons decays at a distance 
from the interaction vertex > 10 m, and is therefore reconstructed in the muon system. 
For the events thus selected, we define N 2 as the number of events in the subsample where 
a second particle with transverse momentum > 10 GeV is identified in the muon system. 
The search for the second particle should be as inclusive as possible, in order to minimise 
the corrections to the ratio 

which is a function of the slepton lifetime. Its dependence on the NLSP lifetime cr in 
meters is shown in Fig. [| for four out of the eight sample models. 
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Figure 3: The ratio R = N^/N\ defined in the text as a function of the slepton lifetime 
ct. The curves corresponding to the model points 1, 5, 6, 8 are shown. 

The probability for a particle of mass m, momentum p and proper lifetime r to travel 
for a distance L before decaying is given by P(L) = e~ mL / pcT . _/y 2 [ s therefore a function 
of the momentum distribution of the slepton, which is determined by the details of the 
SUSY spectrum. One needs therefore to be able to simulate the full SUSY cascade decays 
in order to construct the ct-R relationship. 

Relevant for the precision of the SUSY breaking scale measurement is the error on the 
measured ct. This can be extracted from the curves shown in Fig. 0. 

The precision calculated according to this formula is shown in Fig. ^ (left side), for 
model sample point 1 and an integrated luminosity of 30 fb _1 . The full line in the plot 
is the error on ct considering the statistical error on R only. The available statistics is a 
function of the strongly interacting sparticle mass scale. 

We parameterise the systematic error as a term proportional to R, added in quadrature 
to the statistical error. We choose two values, 1%R and 5%R, and propagate the error to 
the ct measurement. The results are represented by the dashed and dotted lines in Fig. [| 
(left side). 

For the models with squark mass scales up to 1200 GeV, assuming a 1% systematic 
error on the measured ratio, a precision better than 10% on the ct measurement can be 
obtained for lifetimes between 0.5-1 m and 50-80 m. If the systematic uncertainty grows 
up to 5%, the 10% precision can only be achieved in the range 1-10 m. If the mass scale 
goes up to 2 TeV, even considering a pure statistical error only, a 10% precision is not 
achievable. However a 20% precision is possible over ct ranges between 5 and 100 m, 
assuming a 1% systematic error. 

Using the measured values of ct and the NLSP mass, the SUSY breaking scale ^/F can 
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cx(m) VF/100TeV 



Figure 4: Fractional error on the measurements of the slepton lifetime ct (left side) and 
the SUSY breaking scale (right side), for model sample point 1 only. We assume an 
integrated luminosity of 30 fb" 1 . The curves are shown for three different assumptions 
on the fractional systematic error on the R measurement: statistical error only (full line), 
1% systematic error (dashed line), 5% systematic error (dotted line). 

be calculated from Eq. (|j), where B = 1 for our case where the NLSP is a slepton. The 
fractional uncertainty on the VF measurement can be obtained adding in quadrature one 
fourth of the fractional error in ct and five fourths of the fractional error on the slepton 
mass. In Fig.Jl (right side), we show the fractional error on the yr measurement as a 
function of yF for our three different assumptions on the ct error. The uncertainty is 
dominated by ct for the upper part of the yF range, and grows quickly when approaching 
the lower limit on y/F. This is because very few sleptons survive and the statistical 
error on both m| and ct gets very large. The error on y/~F is better than 10% for 
1000 < VF < 4000 TeV. 

4 Conclusions 

We have shown how the ATLAS detector at the LHC can be used to detect and measure 
the mass and lifetime of massive long-lived charged particles produced in the framework 
of supersymmetric models with gauge-mediated SUSY breaking, where a slepton is the 
NLSP and decays into a gravitino with a lifetime in the range 0.5 m < ctnlsp ~ 1 km. 
At the LHC, a large amount of SUSY particles of this kind is generated, while the SM 
background can be reduced at a very low level with appropriate cuts. 

The lifetime measurement can be performed by a counting method with a good accu- 
racy, allowing the computation of the value of the SUSY breaking scale yF. 

The peculiar problems of the MSEP detection, related with the event trigger and the 
time-of-flight measurement, have been discussed. The delicate point of the possibility of 
measuring the specific ionisation by the MDT's and by different sectors of the hadronic 
calorimeter has been addressed. We stressed the importance that the specific ionisation 
measurements have for a clear MSEP identification. 
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